AOCS Performance and Stability Validation for a 
160 -m Solar Sail with Control-Structure Interactions 

Stephanie Thomas* and Michael PaluszeM 
Princeton Satellite Systems, Princeton, NJ 085^2 

Bong Wie* David Murphy** 

Arizona State University, Tempe, AZ 85287 ATK Space Systems, Goleta, CA 93117 


Future solar sail missions, such as NASA’s Solar Polar Imager Vision, will require sails 
with dimensions on the order of 50-500 m. We are examining a square sail design with 
moving mass (trim control mass, TCM) and quadrant rotation primary actuators plus 
pulsed plasma thrusters (PPTs) at the mast tips for backup attitude control. Quadrant 
rotation is achieved via roll stabilizer bars (RSB) at the mast tips. At these sizes, given 
the gossamer nature of the sail supporting structures, flexible modes may be low enough 
to interact with the control system, especially as these actuators are located on the flexible 
structure itself and not on the rigid core. This paper develops a practical analysis of the 
flexible interactions using state-space systems and modal data from finite element models 
of the system. Torsion and bending of the masts during maneuvers could significantly affect 
the function of the actuators while activation of the membrane modes could adversely affect 
the thrust vector direction and magnitude. Analysis of the RSB and TCM dynamics for 
developing high-fidelity simulations is included. 

For control analysis of the flexible system, standard finite-element models of the flexible 
sail body are loaded and the modal data is used to create a modal coordinate state-space 
system. Key parameters include which modes to include, which nodes are of interest for 
force inputs and displacement outputs, connecting nodes through which external forces and 
torques are applied from the flex body to the core, any nominal momentum in the system, 
and any steady rates. The system is linearized about the nominal attitude and rate. The 
state-space plant can then be analyzed with a state-space controller, and Bode, Nyquist, 
step and impulse responses generated. The approach is general for any rigid core with a 
flexible appendage. This paper develops a compensator for a simple two- mass flex system 
and extrapolates the results to the solar sail. A finite element model of the 20 m solar sail 
by ATK Space Systems, recently validated in ground tests, is used to demonstrate the sail 
analysis approach. 


I. Introduction 

The sail in question is a square sail with four gossamer masts and four triangular sail membranes sus- 
pended between them. IVim control masses move along the masts via lanyard lines and stepper motors, 
providing a center-of-mass/center-of- pressure offset which is used for pitch and yaw control. Stabilizer bars, 
which connect the membrane corners to the mast tips, can be simultaneously rotated to create a windmill 
torque for roll control. The coordinate frame of the sail - shown in a sun-synchronous Earth orbit - is given 
in Figure 1 on the following page. 

The dawn-dusk sun-synchronous (DDSS) orbit is a possibility for a LEO sail validation mission. A sail 
for this mission would likely be on the order of 40 m and would require a controller fast enough to handle the 
solar pressure disturbances. Solar Polar Imager is an example science mission at 0.48 AU and 75°inclination, 
requiring a 160 m sail to achieve the science orbit in about 5-6 years. A sample maneuver from this mission 
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Figure 1. Sail Coordinate System 


is a 3 hour, 70° rotation during the inclination cranking phase; this maneuver happens twice per orbit. In 
this orbit a very low control bandwidth of 0.002 rad/s is used. 5 

The software used to develop this paper is part of the Solar Sail Control Toolbox (SSCT), a MATLAB- 
based toolbox for analyzing and simulating sailcraft control systems. 4 The package is built around Princeton 
Satellite Systems standard CAD package and disturbance models, with specialized dynamics models as 
needed. Attitude models include a two-body gimbaled boom model and two moving mass models - one 
with multibody dynamics and one assuming stepper motors are the actuator. Sail models include a flat 
plate, a sphere, and a square sail with multiple actuators. Analysis demos are included which study the 
use of moving masses and roll stabilizer bars in detail. Simulations include near-Earth validation orbits 
and heliocentric missions such as Solar Polar Imager. The modular nature of this toolbox allows users to 
study various aspects of the sailcraft control problem, from orbit maneuvers to attitude control validation. 
The toolbox facilitates both simplified, first-level analyses of control concepts and high-fidelity simulation 
with full ephemeris and gravity models and numerical disturbance analysis including planetary albedo and 
combined thermal/optical models of the solar radiation pressure force. Third-party tools, such as structural 
analysis and trajectory optimization routines, can be incorporated into SSCT-based analyses since all source 
code is provided and simulations are script-based. Since MATLAB scripts are flexible, a single CAD model 
can serve as a database for both PSS simulation data and third-party tool inputs. 


II. Roll Stabilizer Bar Analysis 


Initial analysis of the roll bars assumed that they produce a simple windmill torque. However, high 
fidelity simulations of a full CAD model with disturbances in SSCT have indicated significant coupling of 
the yaw and pitch axes during roll stabilizer bar (RSB) actuation. An analysis of the complete geometry, for 
any sun vector to the sail, reveals the source of the coupling to be related to the projected area differences 
between the sail quadrants, which must be compensated for by the other actuators. 

For the solar sail, the nominal location of the sun vector in the body frame is the +x axis. In general 
there will be azimuth and elevation angles a and /? between the sun vector and +x due to thrust vector 
control commands or pointing errors. The sun unit vector is then expressed as: 


u 


sun — 


cos a cos /? 
sin q cos /? 
sin/? 
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Each sail quadrant can be twisted about its tack line using the RSBs. The geometry is shown in Figure 2. 
The sail quadrants are in the Y/Z plane, with the masts along the Y and Z axes. The four stabilizer bars 
are always rotated the same angle S. 



Figure 2. RSB Geometry 


The centroid of a triangle is the average of the vectors to its corners. Expressed in the body frame, the 
sail centroid vectors are 

[ 0 0 0 0 

-♦ Imast 

T cm — ^ 1—1 — 1 1 

[1 1-1-1 

where each column is the centroid for one quadrant, starting from the +Y j+Z quadrant and moving coun- 
terclockwise. Imast is the length of the mast from the sail center to the tip. We see that for every quadrant, 

Item I = y/2 ^ 


The geometric relationship between the RSB angle 5 and the resulting sail angle 9 is: 


sin - = — 

2 


0 


2 ^ bar sin 2 

I/ 

2 L sail 

0 ^bar • & 

- 2 T ~, s ' n 2 


as shown in Figure 3 on the following page. Zf, ar is the total length of the RSB and l sa ii is the hypotenuse of 
the quadrant. The rotated stabilizer bar and sail each make an isosceles triangle which share an edge. The 
two approximations used are the small angle approximation for # and Z(, ar /2 Imast- For this definition of 
the sail angle, a positive RSB rotation results in a negative sail angle. 

The sail normals are defined in terms of 9 as: 


Tlsail — 


cos# 


cos# 


Tf.sin# 75 sin # 
-72 sinf? 72 sine 


cos# 

-72 sm9 

72 sin0 


cos# 

-^sin# 

“75 sin0 


The sun angle experienced by each sail quadrant is found from the dot product of its normal with the 
sun vector. 


')sun — COS 


= COS 


( "O' sun ’ f^saiz) 


1 


(cos a cos 0 cos 9 H — j= sin 9 
\/2 


( sin a cos /3 — sin /? \ 
sin a cos (3 + sin 0 
— sin a cos 0 + sin /? 

V — sin a cos 8 — sin 0 ) 
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2D trigonometric relationships 



Notational 3D geometry 


Figure 3. Relationship between RSB angle and sail angle 


Thus, the result of twisting the sail quadrants in the presence of nonzero a and @ sun angles is a different 
sun angle for each quadrant. Figure 4 indicates the sun angle each quadrant experiences for a small set of a 
and (3 angles and a full range of RSB angles. 




Figure 4. Example sun angles 


The roll torque generated by rotating the sails, given a total solar pressure force F sun aligned with the 
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X axis, can be simply approximated as 


T = -F sun V 2^ sind = F sun V2 - sin | 

since l m ast/ 3 is the moment arm to the centroid of the sail quadrant. However, for a nonaligned sun vector, 
the torque must be computed from 

F — 'y ) Term x F Si 

i 

The solar force on each sail quadrant (excluding thermal effects) is 

F Si = p s Acosy sun (2(a s cos'f sun + a d /3)h + (a a + a d )u s ) 


where 

O a + + o'd + a t — 1 

Since the force depends on the effective area of each sail, the slightly different angles for each quadrant 
with respect to the sun vector results in a different force magnitude on each sail. Therefore while one might 
expect a pure roll torque from RSB actuation despite the angle of the sun vector, in reality the introduced 
mismatch between the sail forces results in torques along the pitch and yaw axes as well. 

Writing out the cross products and summing over the quadrants gives 

— 4sin#(F 5l + F S2 + F S3 + F Si ) 

-jocose (F S1 + F S2 - F S3 - F Si ) 

75 cos 6(-F Sl + F S2 + F S3 - F Si ) _ 

which, for a purely specularly reflective surface ( o a = a d = 0), reduces to a dependence on the combinations 
of cos 2 7*. 

Substituting the sun angle definitions in Equation 1 into the torque equation results in a lengthy equation 
we will not display, but we can easily consider the limiting cases. When /? = 0, we have Fj = F2 and F3 = F 4 . 
In this case, a change in a leads to a nonzero T y , but T z = 0. Similarly, when a = 0, iq = F4, F2 = F3, (3 
leads to a nonzero T z , but T y = 0. When both a and 0 are nonzero, then T y and T z are both nonzero as 
well. 

An example is taken from 40 m square sail in a 6 pm DDSS orbit, given in Table 1. The angle between 
the sun vector and the orbit normal varies with both the orbit period and time of year. The sail is nominally 
aligned with the LVLH frame to prevent gravity gradient torques, with a -90 degree yaw (z) rotation to bring 
the sail from edge-on the sun to approximately full-on. Figure 5 on the next page indicates the sun angles 
(azimuth, elevation, and total angle) and the RSB torque so that the connection between them is clear. The 
torque is calculated for an RSB angle of 10 degrees (sail angle of 0.19 degrees) and sail tack angles of 0 and 
30 degrees, that is, LVLH yaw angles of -90 and -60 degrees. Note that even for a tack angle of zero degrees, 
the sail experiences about an 11 degree sun angle due to the LVLH alignment. 


Table 1. DDSS Epoch and Orbital Elements 


Year 

Month 

Day 

H 

M 

m 

2010 

3 

15 

16 


0 

a (km) 

i (rad) 

fi (rad) 

w (rad) 

e 

M 

7971 

1.7891 

1.5708 

3.022 

0 

0 



From these, one can clearly see that rotating the stabilizer bars creates pitch and yaw torques in addition 
to roll torques when the sail is not perfectly aligned with the sun. These torques oscillate at orbit rate due to 
the oscillation in the LVLH sun angles around the orbit, but the mean value is approximately zero whereas 
the roll torque is constant. Increasing the sail’s tack angle in this orbit to 30 degrees creates a large offset in 
the sun’s azimuth angle (green lines). This is turn creates an offset in the pitch torque, while only slightly 
changing the amplitude of the yaw torque. The small change in yaw torque with tack angle is because, as 
we mentioned previously, T y is more strongly linked to a than T z , and in this case a changes a lot with the 
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Figure 5. An Example of RSB Torque 

tack angle while (3 remains about the same. In general, the geometry of the sun vector, and not just the 
magnitude of the sun angle or the sail angles, drives the resulting three-axis torque distribution. The roll 
torque does remain dominant, indicating that the RSBs are still an effective actuator if sufficient pitch and 
yaw control is also available. The off-axis torques could be directly accommodated by feedforwarding the 
torque estimates to the pitch/yaw controller, improving the overall peformance. 

This calculation has been compared to the full disturbance model, which includes a complete CAD 
description of the components and multiple disturbance types. The torque from the disturbance model for 
a 30°tacked sail, on the right in Figure 6 on the following page, is comparable to the output for a tacked sail 
from this derivation (green line) which uses only solar radiation pressure (SRP), validating this flat-plate 
analysis approach for this model. 

Lastly, we consider a closed-loop simulation with a PD controller and the full disturbance model in 
Figure 8 on page 8. The commanded maneuver is a 30 degree slew in yaw from a a sun-facing attitude. 
The roll axis inertia is roughly twice the other axis, so we see the effects of Euler coupling. We can see that 
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Figure 7. Calculated SRP Forces 


both the roll and pitch axes become involved but that the attitude stabilizes, since we have enough control 
authority in the TCMs. 


III. Trim Control Mass Analysis 

The equations of motion for a two phase stepper motor in ab coordinates (physical, orthogonal winding 
axes) are 1 

= Z ( v a ~ Ri a + K m u)sin(N r 9)) 
f = lK - Rh - K m ujcos(N r e)) 

% = 7 ( K m i b cos(N r 8 ) - K m i a sin{N r e ) -Bui) U 

dO , 

where i a and it, are the phase currents, v a and v b are the phase voltages which will be used for control, N r 
is the number of rotor teeth, L is the inductance of the windings, R is the winding resistance and K m is the 
motor torque constant in units of Nm/A, u) is the rotor angular rate and 8 is the rotor angle. J is the rotor 
inertia. Each phase winding is assumed to have the same resistance and inductance. The response is shown 
in Figure 9 on the following page. 

The phase voltages switch and the rotor moves to the new position. The stepper motor in effect has its 
own feedback loop to maintain the position so no external feedback position control is required. The stepper 
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motor dynamics are very high frequency so it is not necessary to model the dynamics when analyzing the 
solar sail control system. 

For the purpose of simulations we can assume the trim control masses can attain any speed instanta- 
neously, i.e. they do not accelerate. Angular momentum of the system must be conserved, so when the 
masses are moving the vehicle rates must be adjusted accordingly. An example of a 70 degree pitch ma- 
neuver for a 160 m, 450 kg sail in SPI transfer orbit with 5 kg moving masses 7 is shown in Figure 10. The 
moving mass velocity is limited to 5 cm/s. The disturbance model consists of solar radiation pressure on 
purely specular flat plate quadrants. The controller is a PID with a bandwidth of 0.001 rad/s. 




Figure 10. Moving mass actuation demo for a rigid model in SPI orbit 


IV. Sail FEM Model 

The model under study is the 20 m workhorse sail recently deployed at NASA Glenn Research Center 
Plum Brook Station (Sandusky, OH) by ATK Space Systems (Goleta, CA). 3 A 1028 node finite-element 
model is available for this sail. The first 30 modes range from about 0.5 to 1 Hz. The model and modes can 
be effectively visualized in MATLAB as seen in Figure 11 on the following page. 

One control method under consideration is a suite of pulsed plasma thrusters mounted at the tip of the 
masts. 6 On the +Y axis for this model, the corresponding node at which the force will be applied is 40001. 
A +X force applied here will produce a negative rotation about the Z axis. 

We require a state-space form of the dynamics, i.e. 

x = Ax + Bu 
y = Cx + Du 

Our vehicle consists of a rigid core and a flexible appendage. We begin by defining a set of nodes on the 
appendage. If the vector from system center of mass ( C sys ) to the origin is c, and the nodal displacements 
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Figure 11. Solar sail mode visualization in MATLAB; on the left, mode 1, on the right, mode 22. 
are p fc) then the vectors from C sys to the nodes are d k — c + r k + Pk- This is depicted in Figure 12. 


Flex 



Figure 12. Node vector definitions 

We also define the modal displacements rj, where the transformation matrix $ gives the nodal displace- 
ment per unit modal displacement for each mode. $ is constructed as a 3n by m matrix, where n is the 
number of nodes and m is the number of modes. This can be written compactly in matrix notation as 

P = $V ( 2 ) 

Our goal is to develop state-space equations using the modal coordinates. We begin with the total 
angular momentum of a rigid core with a flexible body about the spacecraft center-of-mass. Note that the 
sum includes the displacement of the core from the system center of mass. 

n 

H = A(Iu) + J2 m h D kDk (3) 

k=0 

where A is the rotation from the body frame to the inertial frame and Dk = Adk is in the inertial frame. 
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Taking the inertial derivative and transforming to the body frame we have 

n 

T = Id) ~iW X Io) ^ ^ m k d k d k (4) 

fc= 0 

where I is the inertia matrix of the core and T is the external torque on the spacecraft. The definition and 
derivatives of d are 

d = r + c + p 
d = u)*d + p + c 

d = u) x u x d-\-di x d-\- 2u ) x (p + c) + p + c 


where c, assuming ro is zero, is 


_ _ y' m k {r k + pk) 

h mT 


Next we write the inertial equations of motion for each nodal mass, including the motion of the spacecraft 
center-of-mass, X, produced by the total external force. We drop the nonlinear terms in u, p and c. 

fk = m k (d k +X ) (6) 

y F 

fk ~ m k (—d k d) + p k + c + ^ — ) (7) 

where f k is the sum of the internal forces exerted by adjacent nodes and any external forces on the nodes. 
The internal mass forces can also be represented as a spring force, giving 

fk = ~Yl k ki ( Pi ~ Pk) + F k (8) 

i 

where k is a physical property of the structure and F k are the external forces. Substituting this force 
definition and the definition of c we have 


m k (p k + c) + Kp = -m k (-d k dj+ — — ) + F k 


m k (p k - y -) + Kp = m k d x d) + 
t— 1 rriT 


m T -m k Fj 

F k — m k y — — 

mr , mr 


.mr-nik .. m iPj \ , xr >x • , m r ~ m k 

f pk - > — ^1 + Kp = m k d k uj - i 

' m.rr> mm ' Tfl r r 


Fk — mk — 

U mT 


m k {ykPk - y^.djPj) + Kp = m k d k L) + v k F k - pk^Fj (12) 

where K is symmetric and full. For three-dimensional displacements, p is 3n x 1 and K is 3n x 3 n. Now we 
assign M to be the coefficients of p and N to be the coefficients of F, and we have in matrix form 

Mp + Kp = Md x d) + NF (13) 

$ for the cantilevered flex modes is normalized so that the modal masses are equal to one and 

4> T M<1> = E 

= fi 2 

Therefore, switching now to modal coordinates and making the above substitutions in Equation 13 we have 

MQjj+KQri = Md x d> + NF 
$ T M$i7 + $ T /r<l>7 7 = ^ r {Md x d) + NF) 

+ ft 2 7? = <f> T {Md x dj + NF) (14) 
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This is comparable to the general form of a modal dynamics equation for a flexible system (including 
damping), 

t) + 2ZQ.f] + fi 2 ?7 = Bu (15) 

where Z is the modal damping matrix and u are the nodal inputs, except that having M will result in the 
free frequencies being higher than the cantilevered frequencies. 

Next we move to development of the torque equation (4). First we write an intermediate step of the 
equation by substituting in the equation for d, neglecting only the nonlinear terms in p and c: 

n 

T = Iuj + uj x Iuj - (ui x uj x dk — d%LJ + pk + c) (16) 

fc= o 


The total inertia It is defined as It = I core — ^2 m kd k d k . We gather terms which contribute to It, note 
that the sum of the coefficients of c is zero and po is also zero, and convert to modal coordinates: 

n n n 

T = [I - +w x [J — (17) 

k = 0 k=0 k = 1 

The torque T is the sum of all external torques, including any forces applied at the nodes. 

T = T c + Y^d£F k (18) 


We can now write both the core and mass equations in matrix form with the terms containing ill and r) on 
the left. 


It Xlfc ^Jc, 1 * ' ' 0fc,m 


UJ 

[J2k m kdk(l>k,l] T 


m 




. Efc m kd k 4>k,m] 


Vm 


— u) x Itui 
- fl 2 r] 


+ 


T 

$ r NF 


The left-hand matrix summation terms can be considered a generalized modal inertia matrix, I mo dai- 
The terms indicate the contribution for each mode, which is a sum over the contributions from each node. 
Grouping these terms into a matrix and explicitly writing out the core and node inputs we have 


It Im 

UJ 


(/ T w) x 0 0 


UJ 


' E 3 d 0 x d x ' 


T c 

/£ 4> T M$ 

rj 

— 

0 -n 2 -2 ZQ 


V 

+ 

0 I>» T 








. V . 


L M ^ J 


NF k 


To obtain a state-space form from this equation, we will have to divide by the inertia matrix on the left 
and add kinematics. In addition, we have left in the ( Iu>) x nonlinear term for now, which will need to be 
linearized. There are two common cases, first the presence of a large bias momentum (such as a momentum 
wheel) and second for a steady body rate. In the case of bias momentum and small body rates, the bias 
drives the nonlinear term and the linearized version is simply 


uj x I ui « —h^ui 


In the case of a steady rate, we linearize around the rate and drop terms of second order. 

uj x Iuj = (u + wo) x I(w + Wo) 

~ (—(IuJo) X + UIq I)bJ + UJq Iujq 


To complete our state-space representation we need to add states for the attitude. We consider three easily 
linearized cases, small Euler angles, small quaternion deviations, and perturbations from a rotating frame. 
Small 3-2-1 Euler angles can be integrated to first order directly from the body rates, giving 

e 

u> 


6= 0 E ] 
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If the spacecraft is undergoing a steady rotation and we designate the nominal rate vector as w 0 , then the 
angular velocity is simplified with the small angle approximation as follows, 

u ss 6 + ( E — d x )u>o 


Rearranging to solve for 9 we have 


1 = [ -w 0 * E 


e 

u> 


— Wo 


For small angle quaternions, the first component (<? 5 ) is constant and the derivative reduces to one-half the 
body rates for the latter three terms ( q T ), 


Qs 


O 

O 

O 


Qs 

Qr _ 


0 0 \ 


Qr 



L z J 


UJ 


V, Simple Non-Collocated Systems 

Using a PPT at the tip of the solar sail with sensing at the core is a classic noncollocated problem. 
However, compensating such a system near the flex modes is not trivial and bears review. We first analyze 
one dimensional and planar systems before moving to the full sail model. 


A. One-dimensional system 


We begin our control analysis by compensating a simple 1 dimensional two-mass, one-spring system with a 
non-collocated sensor and actuator. 


Xi, Fj 


X2, F 2 


■0- 


miXi + c(±i - ± 2 ) + k(x\ - x 2 ) = Pi 
m 2 x 2 + c(x 2 - ± 1 ) + fc(x 2 - Xi) = F 2 


In the frequency domain, neglecting damping for the moment, the system becomes 


( 20 ) 

( 21 ) 


CNCN 

+ 

cs 

CO 

1 

w 2 

Fi /mi 

wl 

s 2 + w\ 

F 2 /m 2 


s 2 (s 2 + (wf + V%)) 


( 22 ) 


where u >\ = k/m\ and uj 2 = k/m 2 ■ From this we see that for the collocated case F 2 to x 2 , at low and high 
frequencies relative to the flex resonance the system acts like a rigid body, 1/s 2 . The lead from the zero(s) is 
cancelled by the lag from the pole(s) so that the phase after the flex frequencies remains 180. For F 2 to xi , 
however, at high frequencies the denominator approximates 1/s 4 , resulting in permanent phase lag. This is 
easily seen in the open loop response for both cases in Figure 13 on the following page. 

The system properties are selected to give a natural frequency near the lowest sail mode, m 1 is 10 kg, 
m 2 is 1 kg, k is 0.3, c is 0.03, u> n is 0.574 and $ = [0.0953 — 0.953]. 

We first apply a forward gain for the desired transient response, k = 0.01. The result is unstable. We 
then add phase lead at the crossover frequency for the desired damping, with a maximum phase shift of 40 
degrees at w max = l -2 . This pushes the crossover frequency up, so we lower the gain to 0.001 and raise w max 
to 2e-2, which gives a stable system. There is now a tradeoff between amount of phase lead and excitation 
of the flex mode. Figure 14 on the next page shows the open loop and step responses for phase leads of 40 
and 65. 

Next we add a notch filter to reduce the gain seen by the flex mode. This should allow us to use the 
higher phase lead (i.e. more damping) with less excitation of the mode. We set the notch frequency equal 
to the resonant frequency (0.574 Hz) with a gain reduction of 20 dB and a notch width of 0.1 Hz. See 
Figure 15 on page 15. 
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Figure 15. Compensator with a notch filter at the flex mode 


Lastly, we need to add an integrator to drive the offset to zero. We will need this in the solar sail controller 
since the pointing during delta-V changes controls the thrust vector; we essentially want to limit the average 
pointing error over a maneuver within some tolerance. There is a trade between the time constant of the 
integrator and the forward gain; an infinitely high gain would result in a negligible offset, and correspondingly 
adding the integrator should allow for a higher gain. We will use a slow integrator of 3600 s, and add it in 
parallel with the lead compensator and filter. 



States for Plant * Gain * LFI 


500 



0 

c 

) 

20 

40 

60 

80 100 

120 

140 

160 

180 

500 

f 



_____ 






0 

( 

> 

20 

40 

60 

80 100 

120 

140 

160 

180 




M _ „ ; _ _ 

< 

) 

20 

40 

60 

80 100 

120 

140 

160 

180 


jS*. 










i 









10 ( 

) 

20 

40 

60 

80 100 
Time (min) 

120 

140 

160 

180 


Figure 16. Compensator with integrator added in parallel 

This same loop shaping procedure can be used with more complex systems such as the solar sail. 

B. Planar system 

Now that we have established the response and control of the simplest flex system, we check our theoretical 
development for n flex nodes with a rotational planar system consisting of a core and two masses, having 
two flex frequencies. The core mass is 100 kg, mj is 2 kg and rq is 2 m, and m 2 is 1 kg and m 2 is 4 m, 
and the flex displacements are transverse to the node locations, as shown in Figure 17 on the next page. 
For a spring stiffness of 0.5, the cantilevered natural frequencies are 0.38268 and 0.92388 rad/sec. If the 
flex body is correctly linked to the core, we expect that the response to a force at r 2 will be 4 times the 
response to a unit torque on the core. We see in Figure 18 on the following page that we do get the collocated 
and noncollocated responses we expect. Note the upward shift of the flex frequencies from the cantilevered 
frequencies to the free frequencies. 
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Time (sec) Time (sec) 

Figure 18. Response for a core with a 2 node flex body 

VI. Solar Sail Analysis 

In the case of the 20 m sail, the flex frequencies are relatively high and a PID controller would be sufficient 
for actuation at 0.002 rad/s bandwidth. 7 For larger sails, the modes will scale as some combination of the 
boom modes, a 1/L 2 relationship, and the membrane modes, a l/L relationship, so that the flex modes and 
bandwidth will be considerably closer. There is also the complication that there are many modes very close 
together, so that with overlapping poles and zeros we may not see the clear noncollocated response in the 
simple models. 

First we check the open- loop sail plant with our desired force input at a boom tip. We choose a node on 
the +Y boom and consider forces in the +X and +Z directions, which will produce respectively — T z and 
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Figure 19. Frequency Response including Modes 1-8 


We are interested in developing a controller which may have structural interactions as may occur for 
larger sails. We apply a PID controller at different control frequencies, again keeping modes 1-8 in the plant, 
and check the eigenvalues for stability. For w n equal to 0.001 and 0.01 rad/sec, the PID controller is stable 
with the noncollocated actuation at the tips; at 0.1 however it is not, so this provides an opportunity to try 
the loop shaping procedure outlined in Section V on page 13. 

We will compensate the roll axis since it has an interesting noncollocated response. We begin with a 
gain of 0.6, max phase frequency of 0.2, and a phase gain of 65 degrees. Then we add a notch filter at the 
apparent flex frequency, about 24 rad/sec, with a half-notch width of 5 rad/sec and a gain drop of -20 dB. 
See Figure 20 on the next page and Figure 21 on the following page. 


VII. Conclusions 

We have analyzed the proposed trim control mass (TCM) and roll stabilizer bar (RSB) actuators and 
concluded that it is not necessary to simulate the motor dynamics of the TCMs. The RSBs can produce the 
needed roll control torques but will require compensation for the off-axis torques induced by the sun vector 
geometry. 

We have demonstrated that a relatively simple controller can be used to control a flexible system, including 
the solar sail with tip-mounted actuators. So far we have used idealized actuation and assumed damping in 
the system. In future work, we will incorporate realistic damping, actuator, and sensor models into models 
of . We would perform simulations with the linearized flex plant developed here, a nonlinear core model, 
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+T X . The inertia about the X axis is twice that of the other axes, and the boom tip is about 14.4 m from 
the core center. This is a complex system, and we do not see a simple noncollocated response in both cases 
when considering the first 8 modes. For F x , we see a rigid body response (Figure 19), but for F z , which 
interacts strongly with only one of the 30 given modes (773), we do see a phase drop. 


Input I50062Z to Output tfiX 




Input: <501 ZSy to Output'. FK Input. 1501 25y to Output: IKL Input: I50125y to Output: IKLF 



Figure 20. Frequency response of compensated sail (including modes 1-8), with a pure gain, gain with lead 
compensation, and compensator with notch filter. 
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Figure 21. Step response of compensated sail (including modes 1-8) 

nonlinear sensor and actuator models to demonstrate large angle slews and other maneuvers. 

In analyzing the inputs at a single node we have shown that the structure flex interactions are complex 
and do not necessarily follow simple models. The sail has many closely-spaced frequencies and deciding 
which modes to include in higher-fidelity studies will require detailed analysis. 
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